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1. Introduction 

CP violation was first establislied in — )• 71+ tt^ decays in 1964 [|l|] and since then it has been 
observed in several K and B meson weak decays 1^]. In the standard model of particle physics (SM) 
CP violation in the quark sector of the weak interactions arises from a single irreducible phase 
in the Cabibbo-Kobayashi-Maskawa (CKM) matrix [||] that describes the mixing of the quarks. 
The unitarity of the CKM matrix imposes a set of relations among its elements, including the 
condition VudV*/^ + VcdV*/^ + VtdVjl = that defines a unitarity triangle in the complex plane, shown 
in Fig. [U Many measurements can be conveniently displayed and compared as constraints on 
sides and angles of this triangle. CP violation is proportional to the area of the unitarity triangle 
and therefore it requires that all sides and angles be different from zero. The angle 7 = 03 = 
^E,{~VudV*i^/VcdV*ij) is at present the most difficult to measure. 

An important goal of flavor physics is to over-constrain the CKM matrix. The reason is 
twofold. First, it is desirable to determine its elements as precisely as possible because their values 
are fundamental parameters of the SM. Second, new physics (NP) effects could manifest them- 
selves as inconsistencies between two or more measurements of the CKM parameters [Q]. The 
angle 7 can be measured in decays mediated by tree amplitudes, such as B ^ D^^K: assuming that 
NP does not change the tree-level processes, its determination is not affected by effects beyond the 
SM' and together with the measurement of \Vuh/Vcb \ provides a constraint that can be compared 
with those potentially sensitive to NP. 
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Figure 1: Graphical representation of the unitarity constraint Vi,t/Vj*^ + VcdV*h + ^tdVtl = as a triangle in 
the complex plane |Q| . 

In the following we will describe and discuss the main techniques to measure the angle 7 based 
on the measurement of B ^ DK decays and we will summarize the results. 

2. Measurement of 7 with flavor-tagged B — )■ DK decays 

The most powerful method used so far exploits the interference between b cus and b 
ucs amplitudes in B ^ DK decays, whose relative weak phase is 7. In charged B decays the 
interference is between the B — )• D^K amplitude and the color- and Cabibbo-suppressed B — )• 
D^K amplitude, when D° and decay to a common final state. The leading interfering diagrams 
are shown in Fig. ^ 

NP may appear in O" - d" mixing but the effect is expected to be small and can be taken into account |§. 
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Figure 2: Main Feynman diagrams contributing to tiie B — s> DK decay. The left decay proceeds through a 
b cus transition, while the right diagram proceeds via a — s> ucs transition and is both color- and Cabibbo- 
suppressed. 

Let us introduce the decay amplitudes A(S" D^K^) =Ace'^^, A{B^ D^K^) =Aue'^^"^'^\ 
A{D° ^ f)= Afe'^f and A(D° ^ /) = Aje'^f, where / is a generic final state of the D meson. 
The parameters 5c, 5u, df and 5j are strong phases, and 7 is the weak phase difference between 
B — )• D^^K and B D^K . A^, Ai„ Af and Ay^ are real and positive. The decay amplitude of 

B ^ [f]DK- is 

A{B- ^ [f]DK-) = A,A/e'(^^+^/) + A„Aje'(^"+^/-^) , (2. 1) 

and the rate can be written as 

r(B- ^ [f]DK-) = aIa] (a) /a] + 4 + 2rBA//Ay-Re(e'(^''+^°-^))) , (2.2) 

where rg = Au/A^, 5b = 5u — 5c and 5d = 5j — 5f. The rate of the charge-conjugate mode is 
obtained from Eq. ^ by replacing 7 with —7. Since with excellent approximation there are no 
other contributions besides the ones in Fig. ^, 7 and the unknown hadronic parameters and 5b can 
be constrained in a theoretically clean way by measuring the yields of B^ and B , provided that the 
D decays are chosen appropriately. The same concept and equations apply to D*^K^, B^ — > 

D^K*^ and to the flavor-tagged — )• D^K*^^ decays with small, though important, modifications!^. 

The parameters and 5b are not expected to have the same value in different B D^*)^K^*) 
decays. The magnitude of is crucial because it measures the interference which allows the 
extraction of 7. For charged B decays rg Ki Cf {VcsV^b/'^usVch I ~ 0. 1 — 0.2, where cp ^ 0.2 — 0.4 is 
a color suppression factor. For neutral B decays rg can be as large as 0.4, but the typical yields are 
significantly smaller making these decays less competitive at present. 

In the next three sections the main experimental methods explored so far are briefly discussed. 

2.1 The GLW method 

In the method proposed by Gronau, London and Wyler (GLW) the D meson is recon- 
structed to CP eigenstate final states fcp±, such as K+K {CP = +1) or K^K° {CP = -1). There- 
fore Ay/Aj = 1, 5/) = 0, 71 for CP = 1,-1, and Eq. becomes 

r{B ^ [fcp±]DK-) =AlAl^Jl+rl±2rBCOsi5B-Y)). (2.3) 

^The decay B — > DK* can also be measured as part of the Dalitz plot analysis of S — > DKn. See |^ and references 
therein. 
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Equation 2.3 and its CP-conjugate are used to define four observables which depend on the physical 



parameters 7, 5b and rg: 

^ _ nB-^D'lp^K-)-T{B+^D^CP^K+) _ ±2rg sin sin y 

T{B ^ D^p^K- ) + r(B+ ^ DIp±K+ ) 1 + r| ± 2rB cos 5b cos 7 ' 

Rcp± = 2riB-^D^K-) = 1 + 4 ± cos 5^ cos 7, (2.5) 

where r{B- D%^K-) = r{B [fcp±]DK-)/BF{D^ fcp±) and r{B- D^R-) = ~ 
T{B- [K-k+]dK-)/BF{D^ -^K-k+). Acp± and Rcp± are bound by the relation Rcp+Acp+ + 
Rcp-^cp- = 0. 

Both BaBar and Belle have reconstructed B D^K and B -> D*^^K^ decays with D*° 
D°7r°. BaBar has also selected the mode D*° D°7 and the decay B D'^K* with K*^ 
K^s^^- mesons have been reconstructed in CP-even (K^K^ and 71+ tt^) and CP-odd (^^tt", 
K^^ and K^co) eigenstates. Recendy also CDF has measured B — )- D'^K with D° ^ K^K 
and 71+ 71^. The data samples used by BaBar and Belle consist of 382 and 275 million BB pairs, 
respectively, while CDF has used 1 .0 fb^ ' of data. The results are summarized in Table |I| 

Due to the discrete 8-fold ambiguity in the extraction of 7 from Rcp± and Acp±, and the still 
large uncertainty on rg, the GLW measurements have a poor constraining power on 7 when they are 
considered alone. However, they generally improve the knowledge of rg, 7 and 5b when combined 



with the results of the Dalitz method discussed in Sec. This aspect is further discussed in 
Sec. I 



Table 1: Summary of Rcp± and Acpi measurements. 



Mode 


Experiment 


AcP+ 


Acp- 


1^CP+ 


Rcp- 




BaBar [JlC 
Belle [In 

CDF ni2 


P 

] 
] 


0.27 ±0.09 ±0.04 
0.06 ±0.14 ±0.05 
0.39±0.17±0.04 


0.09 ±0.09 ±0.02 
0.12±0.14±0.05 


1.06±0.10±0.05 
1.13±0.16±0.08 
1.30±0.24±0.12 


1.03±0.10±0.05 
1.17±0.14±0.14 


B D*^'K 


BaBar ^] 
Belle [111 


0.11±0.09±0.01 
0.20 ±0.22 ±0.04 


0.06 ±0.10 ±0.02 
0.13±0.30±0.08 


1.31±0.13±0.03 
1.41±0.25±0.06 


1.09 ±0.12 ±0.04 
1.15±0.31±0.12 


B o*'a:* 


BaBar [ U 


] 


0.09 ±0.13 ±0.06 


-0.23±0.21±0.07 


2.17±0.35±0.09 


1.03±0.27±0.13 



2.2 The ADS method 



In the method proposed by Atwood, Dunietz and Soni (ADS) 015| ] the D meson is recon- 
structed in doubly Cabibbo-suppressed (DCS) states. We consider D" K^n^ as an example in 
the following discussion. The decay rate of the process B — [K^71^]dK^ is the result of the inter- 
ference between B D^K^ followed by the DCS D" K^k^, and the suppressed B D^K 
followed by the Cabibbo-allowed D° — )• K^n^. From Eq. 22 we find 



r(g+ ^ [^+7r+]z)^+) 



rl + rl + IrBVD cos(5b + 5/) T 7) > 



(2.6) 
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where both ro=A f/Aj = |A(D° K^K ) /A(D° K 71+) | and the strong phase difference 5d 
are measured independently |§, [T6|]. Defining Rads and Aads as 

r(g- ^ [/:+7r-]Djr-) +r(g+ ^ [/:-7r+]o/:+) 



r{B [K-7i+]DK-)+r{B+ [K+k-]dK+) ' 
r{B [K+7i-]DK-)-r{B+ [K-n+]DK+) 
r{B [K+7i-]DK-)+r{B+ [K-k+]dK+) 



(2.8) 



it follows 



Rads = rl + rj) + Irs r^cos /cos (5b + 5d), (2.9) 
Aads = 2rB r/j sin 7sin(5s + 5d)/Rads ■ (2- 10) 



Since ro^KK) = (5.78 ±0.08)% ^\J\\ and rg is expected to be around 10%, the interference effect 
can be quite large. Similar relations, with small modifications, are derived for B D*^K and 
B~ D^K* decays (see ^ and refs. in Table |) and for multi-body D° final states ^ Both 



BaBar and Belle have reconstructed the decay B D^K followed by D° ^ /T+tt on datasets 
of 467 and 657 million BB pairs, respectively. BaBar has also selected B — )■ D*^K with D*" — )• 
DOtt" and D*° D^y (467 x 10*^55), and B D^K* with K* K^K- (379 x 10^ BB). On a 
dataset of 465 million BB pairs BaBar has performed the first measurement of flavor-tagged decays 
B° D^K*^ with K*° K+k', selected in the D" final states K+K-, K+n-K^^ and /:+7r"7r"7r+. 
The results are summarized in Table ^. Due to the smallness of the involved branching fractions no 
evidence of signal has been observed so far and the null measurements have been used to set upper 
limits on rg. The strongest hint of signal has been reported by BaBar with a statistical significance 
of 2.6a for B+ — )• D^[K^n'^]K^ . Even when a signal is observed the constraining power on /of the 
ADS method is weak when it is used alone, but it becomes significant when the ADS information 
is combined with the other methods. This aspect is discussed in Sec. ||. 

2.3 The Dalitz or GGSZ method 



If decays to a 3-body final state such as Z) — K^tl^il^, the decay amphtudes of and 
can be written d&Afe^^f = f{m^_^m\) and Ay^e'^'' = f{m\,m^_), where \ 
masses of the A'^'tt^ and K^n^ combinations. The rate in Eq. 7J1 becomes 



T{B^ ^[KIk-K+]dK^)o. \f{ml,ml)\^ + rl\f{ml,ml)\^ + (2.11) 

2rB |/(m^ , m| ) I |/(m| , m^) I cos (5b + 5d {m% , m| ) =F 7) , 

where 5D{m^,m\) is the strong phase difference between f{m\,m^) and f{m^,m\). The ampli- 
tudes f{m\,mL) are measured through a Dalitz plot analysis on a large sample of flavor-tagged 



decays. Therefore, the B yields in Eq. 2.1 1 only depend on the unknowns y, Sb and rg. The great 
advantage of this method [^] is that the relation between the signal yields and the physics param- 
eters varies over the Dalitz plot, making possible the extraction of 7 with only a 2-fold ambiguity 
(7 — 7> 7+ 180°). Furthermore, the K^n'^K^ branching fractions is relatively large (~ 3%). 



Since the direct extraction of rg, 5b and 7 through a maximum likelihood fit (MLF) using Eq. [2.11 
overestimates rs and underestimates the statistical error of 7 and 5b, it is convenient to express 
Eq. [2.11| in terms of the cartesian coordinates x± = Re[rBe'^^''^^^], y± = Im[rBe''^'*^^''], 

r(BT ^ £,0j^ kIk+k-]K^) oc + (4 +4)|/i|2 + 2 [x^Re[/^/i] +3;^Im[/^/;^]] (2.12) 
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Table 2: Summary of Rads A^os measurements, and limits on rs- 



Mode 


Experiment 




^ADS 






BaBar §0 
Belle [|l| 
BaBar [P 


] 


1.36±0.55±0.27 
0.8 ±0.6+2? 

— \J.D 

1.2±1.2±0.9 


-0.70±0.35l[!;'i'^ 
-0.13+;[^I±0.26 

— U.OO 


[0.09,0.193] @95%CL 
<0.19 @90%CL 
<0.19 @95%CL 


B- ^D*^^[D^^7l^^]K- 
B- ^D*°[D°y]K' 


BaBar Q 
BaBar @ 


1.76 ±0.93 ±0.42 
1.3±1.4±0.5 


0.77±0.35±0.12 
0.36±0.94+°;25 


[0.007,0.176] @ 95% CL 


B- -^D'^K*- 


BaBar Q 


6.6±3.1±1.0 


-0.34 ±0.43 ±0.16 


[0.17,0.43] @ 95% CL 
(combined with GLW) 


^K+K- 
^ K+ It- It- K+ 


BaBar |§2 
BaBar |p2 
BaBar |g2 


1 
1 
1 


6.7+^;° ±1.8 
6.0+^;^ ±0.9 
13.7+^'5^±2.2 




[0.07,0.41] @ 95% CL 



where the notation was simphfied using f± = f{m^,m^). The extraction of x± and y± with a MLF 
using Eq. 2. 12 is unbiased. The physics parameters 7, and 5^ are extracted from x± and y± with 
a frequentist statistical procedure. 

BaBar has used this approach to measure the angle 7 with B — ^ D^^K , D*^K (D**' — ^ D°7r" 
and D°7) and D'^K* {K*- K^K^), with D° K^K+K- and K^K+K- using a sample of 386 
milUon BB pairsQ. Belle has selected B D**)"/:" decays (D*° D^K^ and D°7) with D° 
K^K+K- on a sample of 657 milhon BB pairs [||]. The D° ^ KIk+k- and D° ^ K^K+R- 
decay amplitudes are determined with Dalitz plot analyses of large and very pure samples of D" 
mesons from D*+ — )• D°7r+ decays produced in e+e^ — >• cc events. The amplitudes are described 
using an isobar model, consisting of a coherent sum of two-body amplitudes (parameterized using 
relativistic Breit-Wigner) plus a "nonresonant" term. BaBar has described the nn and Kn S-wave 
amplitudes in D° — > K^ti^k^ using a A'-matrix formalism (see [^J] for a detailed discussion). 

The results for x± and y± are reported in Table ^ From the [x±,y±) confidence regions 
both BaBar and Belle determine the la confidence intervals of 7, Sb and rg using a frequentist 
procedure. BaBar finds 7 = (76 ± 22 ± 5 ± 5)° and Belle 7 = (76+{^ ± 4 ± 9)° (B^ ^ D^*^K^), 
where the solution closest to the SM average has been quoted. All results are reported in Table ^. 
Figure ^ shows the {x±,y±) contours for B — > D^K as measured by Belle (top-left) and BaBar 
(bottom-left), the projections of confidence regions onto the (7, r^) and (7, 5b) plane obtained by 
Belle, and the confidence-level as a function of and 7 found by BaBar. The combined B — > 
£)(*)0^(*) measurements of BaBar and Belle correspond to 3.0a and 3.5a evidence of CP violation, 
respectively. 

It is interesting to ask why the statistical error measured by Belle is about twice smaller than 
what BaBar finds, despite the fact that the experimental observables x±,y± have similar uncertain- 
ties. The error on 7 scales roughly as l/rg. Since Belle measurements have central values of 
between 1.5 and 3.5 times larger than BaBar values, though consistent within the errors, the result- 
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Table 3: Summary of x± and y± measurements. The third error is the systematic uncertainty associated to 
the DaUtz model of the D final state. 



Mode 


Experiment 


x+ [10-^] 


y+ [10-^] 


x_ [10-2] 


y- [10-2] 


D — 7 U A 

Dalitz 
GLW 


BaBar 
Belle ||f 
BaBar 


IC 




-6.7±4.3±1.4±1.1 
-10.7±4.3±1.1±5.5 

-9±5±2 


-1.5±5.5±0.6±0.8 
-6.7±5.9±1.8±6.3 


9.0±4.3±1.5±1.1 
10.5±4.7±1.1±6.4 
10±5±3 


5.3±5.6±0.7±1.5 
17.7±6.0±1.8±5.4 


Dalitz 
GLW 


BaBar 
Belle |||, 
BaBar 




13.7±6.8±1.4±0.5 
13.3±8.3±1.8±8.1 
11±6±2 


8.0±10.2±1.0±1.2 
13.0±12.0±2.2±6.3 


-11.1±6.9±1.4±0.4 
2.4±14.0±1.8±9.0 
0±6±2 


-5.1±8.0±0.9±1.0 
-24.3±13.7±2.2±4.9 


B D"K* 
Dalitz 

GLW 


BaBar 
Belle 1 
BaBar 


i 

28 

F 


1 
1 


11.3±10.7±2.8±1.8 
-W.5t\lj±0.6±&.8 
18±14±5 


12.5±13.9±5.1±1.0 
-0.4+]^;^ ±1.3 ±9.5 


11.5±13.8±3.9±1.4 
-78.4+^^-^±2.9±9.7 
38±14±5 


22.6±14.2±5.8±1.1 
-28.ll33;5±4.6±8.6 



ing 7 uncertainty of Belle is smaller. The concept is illustrated in Fig. |^. Fluctuactions of the 7 
error due to the "l/rs" effect are expected to decrease as the relative uncertainty on rg will become 
smaller. 

At present the Dahtz method has the best sensitivity to 7. However, the uncertainty associated 
to the Dalitz model of the D final state is already the dominant contribution to the systematic error 
and it may be difficult to reduce it greatly in the future. To bypass this limit a model-independent 
analysis is required [^]. At the price of a small loss of statistical power [p9|], the method is free of 
model-dependent assumptions on the D decay and therefore it is a promising approach to follow at 
LHCb and at the next generation B-factories. It requires the use of entangled *P(3770) — > DD de- 
cays at tau-charm factories such as CLEO-c, BES-III or next generation Super Flavor factories [Q]. 
It has been recently shown that using 818 pb^' of CLEO-c data an error on 7 of ~ 2° associated 
to the knowledge of the relative strong phase difference of ^ K^n^n^ and ^ K^n^n^ can 
be obtained pO|]. 

BaBar has applied the Dalitz method also to flavor-tagged D^'iT*" decays, with K*^ 
K^n^ and — )• K^n^n^ . On a dataset of 371 million BB pairs 39 ±9 signal candidates have 
been selected. Using the D" Dalitz plot model obtained in the charged B analysis and imposing 
an external measurement of r^, a loose constraint on 7 was set [|T|]. Self tagging B^ — )• D^K*^ 

Table 4: Measurement of 7, 8b and in B — s> Z)^*^*'/^'*'- decays reconstructed by BaBar and Belle. 



Parameter 


BaBar 


Belle 


7 

rB{D'K) 
Sb{D'K) 
rB{D*'^K) 
8b{D*''K) 

TBiD'^K*) 

5b{D''K*) 


76±22±5±5 
0.086 ± 0.035 ± 0.010 ± 0.011 

(1091^^ ±4 ±7)° 
0. 135 ± 0.05 1 ± 0.011 ± 0.005 

(2971^8 ±5 ±4)° 

0.163^1^;°^^ ±0.037 ±0.021 
(1041^3 ±17 ±5)° 


76l}^±4±9 
0-161^o:o38 ±0.011 ±0.049 
(137.4+}^;° ±4.0±22.9)° 

19fi+° °^2,nni 2+0.062 
^•^^0-0.069 ='=^•^^^-0.012 

(341.9+!^:^ ±3.0±22.9)° 
0.564;°;f^^± 0.041 ±0.084 
(242.6+2°;2± 2.5 ±49.3)° 
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decays and time-dependent measurements of B° — )• DK^ decays are expected to be powerful tools to 
measure 7 at LHCb and at the next generation B-factories, where their production will be abundant 
and the large interference (r^ ~ 0.4) can be fully exploited. 
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Figure 3: Top row: Belle; bottom row: BaBar. Ic7 and 2a contours in the {x±,y±) plane of — > D'^K^ 
(top-left); projections of confidence regions for the B — > DK mode onto the (/"b,/) plane (top-center) and 
onto the {^B^Y) plane (top-right); 1- and 2-standard-deviation regions in the {x±^y±) plane of B^ — s- D'^K^ 
(bottom-left); 1-CL as a function of rs (bottom-center) and y (bottom-right) for B D^^K, D*^K and D^K*. 



3. Note on the cartesian coordinates 

The output of the GLW and ADS methods can be expressed in terms of the same cartesian 
coordinates x± = rs cos(5e ± 7) and y± = sin(5e ± 7) measured with the Dalitz method. This alter- 
native way to quote the results can be useful when different methods are compared or combined. In 
this section we will derive the alternative parameterization and will discuss some of the advantages 
with respect to the classic observables. 

From the definitions in Sec. we can write Y{B D^p^iT" ) /r(B" D'-^K^) = \+xi + 
y\ lb 2x_ , from which it follows 



(3.1) 



■■x\^y\. (3.2) 
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Figure 4: a): Geometric definition of 7 in the {x±^y±) plane and pictorial description of how its uncertainty 
changes as a function of rg- Circles represent x±,y± error regions for Z?+ (A) and B (▼) and the dashed 
lines delimit the range of variation of ly. b,c,d): 1- and 2-standard deviation regions in the {x±,y±) plane 
using the GLW (b), ADS (c) and Dalitz (d) methods, on a dataset of lab^'. The assumptions used are 
discussed in the text. 



The same relations hold for x+ and x\ +y\ after replacing B with B+. Equation 3J_ has been 
used by BaBar to measure x± with the GLW analysis'*, whose results are reported in Table ^. It is 
interesting that the measurements of x± performed with the GLW and Dalitz methods have about 



the same uncertainty on datasets with similar size. The constraint given by Eq. 3.2 is much looser 
because of the quadratic dependence on x,y and the fact that ^ 1. Therefore, the GLW method 
measures x± quite precisely but not y±, and this is the reason why it can hardly constrain 7 when 
it is considered alone. When combined with the Dalitz method, however, the overall error of 7 can 
improve significantly. 



The relation was expressed in terms of Rcp± and Aqpj^. 
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Proceeding in a similar way, Eq. of the ADS method can be written as 



that represents two circles in the {x±,y±) plane centered at (— r/jcos^o, r^sinS/)) and with radii 
7?^ = ^ r(gTl) [^T ) "■ It is not possible to determine 7 with only the ADS analysis of B — > 
D^K because the true x±,y± points are delocalized over two circles^. However, the measurement 
can be combined with the Dalitz and GLW analyses to reduce the overall error of x±,y±, and 
therefore the uncertainty of 7. Figure ^ shows the constraints provided by the GLW, ADS and 
Dalitz measurements with B D^K decays on a dataset of lab^^ assuming {y,rB,8B,rD-,8D} = 
{75°, 0.1, 1 10°, 0.06, 191°} in a scenario where the measured observables are centered to their true 
value. The uncertainty on the values of and 5d was neglected when drawing the ADS constraint. 

4. Conclusions 

Despite the fact that BaBar and Belle have collected almost 1.5 ab^' of data, a precise mea- 
surement of the CKM phase 7 is not yet available. This is not surprising if one thinks that the 
typical branching fractions of the involved decays are of the order of 10^^ or smaller and the inter- 
ference term is found to be around 10% in the main decay modes. In fact, the sensitivity reached 
by the B-factories is much better than it was initially foreseen. The decays — > D^*'*'^'^*'^ are 
currently the most sensitive tool. The method where the D meson decays to K^h^h^ {h = 7i,K) 
has the smallest uncertainty, ranging between 15° and 24° including systematics. The dominant 
source of systematic uncertainty currently comes from the Dalitz model of the D final state. 

Using a Bayesian statistical procedure and combining BaBar and Belle results of GLW, ADS 



and Dahtz methods, the UTfit collaboration finds 7 = (78 ± 12)° while CKMfitter using a 



frequentist approach quotes 7= (73^25)° t^l- results are reported in Fig. ^. To reduce the 
error to a few degrees we have to wait for LHCb (a^ ?a 2 - 3° with 10 fb"' [||]) or a Super Flavor 
factory (a^ ^ 1° with ISah-'^ [|]). 
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